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PREFACE 


The  study  reported  herein  was  conducted  using  funds  provided  by 
the  Operations  Division,  Office,  Chief  of  Engineers  (OCE),  under  the 
auspices  of  the  Investigation  of  Operations  and  Maintenance  Techniques 
(IObrC)  program.  Mr.  Milt  Millard  was  the  OCE  Technical  Monitor  for  this 
work. 

The  study  was  conducted  during  the  period  August  1973  to  April  1975 
in  tne  Hydraulics  Laboratory  of  the  U.  S.  Army  Engineer  Waterways  Ex- 
periment Station  (WES)  under  the  general  supervision  of  Messrs.  H.  3. 
Simmons,  Chie’’  o^  the  Hyaroiiics  Laboratory;  F.  A.  Herrmann,  Jr.,  Assis- 
tant Chief  of  the  hydraulics  Laboratory;  Dr.  R.  W.  Whalin,  Chief  of  the 
Wave  Dynamics  Division;  and  Mr.  D.  D.  Davidson,  Chief  cf  +he  Wave  Re- 
search Branch.  The  investigation  was  conducted  and  this,  report  preparea 
by  Mr.  A.  W.  Garcia,  Research  Oceanographer,  and  CPT  F.  C.  Perry,  CE. 

Directors  of  WES  during  this  study  and  the  preparation  and  publica- 
tion of  this  report  were  COL  G.  h.  Hilt,  CE,  and  COL  John  L.  Cannon,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (SI)  units  as  follows: 


Multiply 

fCrct 

miles  (U.  5.  statute) 
cubic  yarns 

pounds  (mass)  per  cubic  foot 
square  feet  per  second 
degrees  (angle) 


§L 

0.30U8 

1.6093W» 
o .76U55U9 
16.0185^6 
0.0929030U 
0.017^5329 


To  Obtain 

metres 
kilometres 
cubic  metres 

kilograms  per  cubic  metre 
square  metres  per  second 
radians 
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BEACH  NOURISHMENT  TECHNIQUES 


A MEANS  OF  PREDICTING  LITTORAL  SEDIMENT  TRANSPORT 
SEAWARD  OF  THE  BREAKER  ZONE 

PART  I:  PURPOSE  AND  SCOPE 

1.  One  raeans  being  considered  for  use  in  the  artificial  nourish- 
ment of  beaches  is  the  artificial  deposition  of  sand  seaward  of  the  surf 
zone.  Saucers  of  this  method  requires  that  the  deposited  material  be 
placed  so  ti:at  it  Is  entrained  into  the  offsho**e  littoral  transport  sys- 
tem with  favorable  probability  of  eventually  being  moved  into  the  surf 

1 2 

zone.  Previous  attempts  to  accomplish  this  * placet,  the  sediment  on 

the  basis  of  grain-size  distribution  and  the  expected  orbita’  velocities 

3 

of  the  wave  regime.  A later  evaluation  of  the  aforementioned  study 
showed  that  there  was  no  appreciable  movement  of  the  deposits,  material 
over  a period  of  four  years. 

2.  Ir  the  referenced  reports,  the  "Sijulstrom  diagram"  was  the 
basis  for  determining  the  boundary  of  eroding  versus  nonerodir.g  veloc- 
ities for  a given  sediment  size.  Although  the  Hjulstrom  diagram  is  ;se- 
ful,  it  apparent.lv  does  not  readily  lend  itself  to  a viable  solution  for 
the  problem  formed  t^d. 

3.  The  objev  ve  of  the  present  stuoy  is  to  develop  a mathematical 
model  for  cotaining  a quantitative  estimate  the  amount  of  sediment 
transported  by  longshore  currents  and  a qualitatj  e estimate  of  the 
amount  of  r.ediment  transported  into  the  sur ' zone.  The  physical  pro- 
cesses go-erning  sediment  mover. nt  within  the  surf  zone  are  not  well 
understood  which  necessitates  abandoning  the  search  for  a quantitative 
estimate  of  sediment  transport  into  the  surf  zone  without  additional 
laboratory  and  field  experiments.  However,  in  an  approximate  fashion, 
some  conclusions  can  be  reached  regaling  the  effects  of  sediment  trans- 
ported into  the  surf  zone.  The  most  favorable  situation  (from  the  view- 
point of  beach  nourishment)  would  be  that  all  artificially  placed  sedi- 
ment move  through  the  surf  zone  and  deposit  on  the  beach.  In  practice, 
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only  & fraction  of  the  artificially  deposited  material  will  be  trans- 
ported into  the  surf  zone.  Eyen  if  this  quantity  does  not  directly 
nourish  the  beach,  at  the  very  least  it  will  help  saturate  the  surf  zone 
with  sediment  and  aid  (at  least  temporarily'  in  retarding  erosion  of  the 
beach. 

b.  From  the  inception  of  this  .study  it  vas  recognized  that  the 
wave  data  would  probably  be  the  least  reliable  input  to  the  model,  and  a 
reliable  verification  of  the  model  would  probably  require  a comprehen- 
sive  field  testing  program.  This  was  pointed  out  as  early  as  1952  by 
several  investigators  who  concluded  that  aiy  study  of  the  movement  of 
sand  by  waves  should  be  accomplished  by  direct  observations,  although  at 
that  tir-.c  the  means  of  doing  so  had  not  bee  . sufficiently  developed. 
Sediment  tracing  techniques  have  advanced  considerably  since  1952,  and  a 
field  verification  of  the  model  described  in  this  report  and  other 
models  under  development  is  being  planned.  In  order  to  provide  at  least 
a tentative  order  of  magnitude  verification,  a wave  hindcast  was  used  in 
the  model  to  compute  the  amount  of  littoral  transport  at  Point  Peder- 
nales,  California.  The  result  of  this  comparison  is  discussed  in  Part  V. 


PART  II:  SUMMARY  OF  PREVIOUS  INVESTIGATIONS 

Description  of  Sediaent  Movement  in  the  Offshore  Zone 

5.  Most  sediaent  transport  produced  by  gravity  vaves  takes  place  ] 

f 

vithin  the  surf  zone.  Turbulence  produced  by  breaking  waves  is  capable 
of  placing  in  action  all  bottom  aaterial  except  for  the  largest  boulders. 

It  is  in  the  surf  zone  that  the  longshore  current  produced  by  vaves 
obliquely  attacking  the  shore  is  most  intense.  Strong  currents  coabined 
with  extreme  agitation  can  carry  a vast  quantity  of  sediaent. 

6.  Seaward  of  the  breaker  line,  conditions  are  not  as  violent, 

and  the  ability  of  gravity  vaves  to  transport  sediaent  is  considerably 

reduced.  In  this  region,  wave  notion  can  be  adtcuately  described  by 
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classical  btm  1 1-aaplitude  wave  theory,  and  the  resulting  sediaent 
transport  can  be  exaained  in  sone  detail. 

?.  Surface  gravity  vaves  start  to  be  affected  by  or  "feel”  botton 
when  they  enter  a depth  approxinately  equal  to  one-half  wavelength.  How- 
ever, it  is  not  until  they  propagate  into  nuch  shallower  depths  that  ap- 
preciable sediaent  novecent  is  produced.  The  depth  at  which  this  occurs 
is  difficult  to  deteraine  due  to  the  large  muaber  of  variables  involved. 

Froc  laboratory  experinents,  the  relative  effects  of  the  various  param- 
eters on  the  initiation,  voluae,  and  direction  of  sediaent  transport  in 
the  offshore  zone  have  been  outlined.^  ’iuaerous  field  observations  in- 
dicate depths  at  which  sediaent  novenenl  is  initiated,  but  all  are  depen-  • 

dent  on  local  conditions.  Occasional  sand  raovment  at  depths  of  100  ft* 
off  La  Jolla,  California,  is  reported  ir.  Reference  7.  Observations  at 
Point  Conception  and  Arguello,  California,  revealed  that  little  or  no 

bottoo  sediaent  movement  occurred  at  depths  greater  than  60  ft  and  only 

8 

occasional  novement  in  depths  between  30  and  60  ft.  Very  precise  aea- 

sureaents  of  changes  in  sand  levels  on  the  shelf  at  La  Jolla,  California, 

q 

revealed  that  depth  changes  of  0.15  ft  occurred  at  a depth  of  70  ft. 

The  occurrence  of  ripples  generated  by  oscillatory  wave  action  at  depths 

* A table  of  factors  for  converting  U.  S.  custoaary  units  of  neasure- 
aent  to  netric  (SI)  units  is  presented  on  page  3. 
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of  150  ft  off  La  Jolla,  California,  was  observed  by  scuba  divers  and 
reported  by  Inman. ^ 

3.  The  node  of  sediment  movement  in  the  offshore  zorc  is  consid- 
erably different  fron  the  surf  zone.  Early  work  by  the  Beach  Erosion 
3oard^  (now  the  Coastal  Engineering  Research  Center)  shows  that  only  in 
the  surf  zone  is  there  an  appreciable  amount  of  sediment  movement  by  sus 
pension.  In  the  offshore  zone,  once  novenent  is  initiated,  sand  moves 
almost  entirely  as  bed  load. 

9.  The  direction  of  sediment  transport  in  the  offshore  zone  is 
considered  to  be  generally  landward.  In  98  percent  of  the  laboratory  ex 
periments  in  Reference  6,  transport  was  observed  in  the  landward  direc- 
tion. This  coincides  with  results  predicted  by  mass  transport  equations 

12 

developed  by  Longuet -Higgins . In  general,  these  equations  predict 

that  under  the  influence  of  gravity  waves  there  are  thin  'ayers  of  water 
which  move  quickly  landward  at  the  surface  and  bottom  wi  h a slowed  sea- 
ward flow  of  water  at  intermediate  depths.  Since  sedime- r transport  in 
the  offshore  zone  is  almost  all  bed  load,  its  corresponding  movement 
will  be  in  the  shoreward  direction. 

10.  Even  though  there  is  a general  landward  movement  of  sediment 
in  the  offshore  zone,  this  movement  is  significantly  affected  by  thi 

longshore  w.u*rent.  Although  describe:,  ns  being  primarily  confined  to 

13  '•*•*  15 

the  surf  zone,  the  works  of  Longuet-Higgins  and  Steward*  and 

Thornton1^  show  that  a significant  longshore  current  can  exist  in  the 

offshore  zone.  Thorntor.'c  verified  this  by  actual  field  measurements  of 

longshore  velocities  and  sediment  transport  ir.  the  offshore  zone. 

11.  In  general,  seuiiw.v.  movement  by  gravity  waves  in  the  off- 
shore zone  can  be  described  as  being  initiated  in  relatively  snallcv 
water  and  moving  primarily  as  bed  load  diagonally  to  the  shoreline  under 
the  influence  of  both  mass  transport  and  longshore  currents. 

Description  of  Theories  on  the  Movement  of  Sediment 
:ioraal  to  the  Shoreline 


12.  Present  theories  on  the  notion  of  sediment  normal  to  the 
shoreline  have  evolved  from  the  null-point  theory,  first  proposed  by 
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Cornoglia,  which  states  that  there  are  two  opposing  forces  being 
applied  to  sediment  on  the  bottom  in  the  offshore  zone.  These  are  the 
onshore  wave  drift  caused  by  unclosed  water  particle  orbits  generated  by 
progressive  waves  and  the  offshore  component  of  the  gravitational  force 
produced  by  bottom  slope.  It  is  theorized  that  a point  of  equilibrium 
exists  between  these  forces  separating  onshore  and  offshore  notion  for 
every  size,  shape,  and  density  of  sedinent  particle.  The  larger  or 
denser  the  sedi lent  particle,  the  closer  to  shore  will  be  its  null  point. 
Shoreward  of  the  null  point  the  par  >icle  will  move  onshore , while  sea- 
ward from  the  null  point  the  particle  will  nove  offshore  until  it 
reaches  a point  where  no  motion  occurs  at  all. 

13.  .Many  recent  investigators  have  concentrated  on  the  null-point 

theory.  Field  studies  at  La  Jolla,  California,  revealed  that  sedinent 

xS 

tends  to  decrease  in  size  in  the  offshore  direction.  This  suggests 

that  sedinent  is  being  sorted  by  the  intensity  of  the  forces  applied  to 

it  and  is  allowed  to  establish  its  point  of  equilibrium  or  null  point. 

19 

Laboratory  experiments  by  Ippen  and  Eagleson  and  Eagleson,  Glenne,  and 
20 

Dracup  indicate  that  although  a null  point  does  exist  for  all  sediment 

sizes  under  all  wave  conditions,  the  net  sediment  transport  at  the  bed 

°1 

is  always  in  the  shoreward  direction.  Johnson  and  Eagleson"  have 
developed  a means  of  estimating  the  onshore  bed-load  transport  by  using 
these  experimental  results.  In  application  of  these  methods  to  a sec- 
tion of  beach  in  southern  California  they  were  able  to  predict  an  aver- 
age annual  onshore  transport  rate  which  at  least  appears  reasonable, 
although  it  cannot  be  verified. 

lU.  Ocher  investigators,  however,  believe  that  sediment  transport 

in  the  offshore  zone  is  entirely  in  the  onshore  direction  due  to  the 

inequality  ol  onshore  versus  offshore  wave-induced  orbital  velocities  in 

13 

the  relation  of  wave  advance  over  a shoaling  bottom.  Inman  indicates 
that  the  onshore  mass  transport  of  water  by  waves  is  continuous  and  is 
offset  by  the  occurrence  of  rip  currents  that  carry  the  buildup  of  water 
in  the  surf  zone  offshore  again.  The  sediment  transport  corresponding 
to  such  a circulation  system  would  have  sediment  moving  onshore  as  bed 
load  under  the  influence  of  waves  and  carried  offshore  as  suspended  load 
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in  a rip  current  where  it  would  he  sorted  normal  to  the  shore  as  the  rip 

op  23 

current  loses  intensity.  Field  studies""*  discount  the  effects  of  the 
rip  current  as  a return  mechanism.  These  studies  indicate  that  there  is 
general  shoreward  movement  of  sediment  as  hed  load  along  the  bottom  with 
the  finer  material  being  driven  into  suspension  in  the  breaker  zone  and 

carried  offshore  at  middepths.  This  also  is  suggested  experimentally 

19  2k 

by  Ippen  and  Eagleson  and  Russell  and  Osorio. 

12.  Other  investigations  indicate  that  the  transport  of  suspended 

sediment  may  be  more  important  in  the  offshore  zone  than  previously  sus- 

25-31 

pected,  particularly  in  the  process  of  sediment  sorting.  However, 

general  net  transport  can  still  be  assumed  to  be  in  the  shoreward  direc- 
tion as  bed  load,  despite  the  increased  evidence  of  offshore  movement  of 
suspended  sediment. 


Description  of  Theories  on  the  Movement  of  Sediment 
Parallel  to  the  Shoreline 

op 

l6.  Galvin0  gives  a comprehensive  review  of  the  theory  and  data 
available  in  North  America  up  to  196?  on  longshore  current  velocities. 
Nine  papers  are  discussed  containing  twelve  equations  of  which  six  can 
be  evaluated  using  available  data.  The  difficulty  in  applying  the  equa- 
tions discussed  by  Galvin  rests  with  the  fact  that  they  give  a mean 
longshore  velocity  and  one  of  the  necessary  requirements  for  accomplish- 
ing the  objective  of  this  study  is  a means  of  calculating  a longshore 
velocity  field.  This  requirement  results  from  the  need  to  know  the 
depth  at  which  a significant  amount  of  the  deposited  material  can  rea- 
sonably be  expected  to  be  entrained  *n to  the  littoral  current.  Gole 
et  al.  (personal  communication)  compared  estimates  of  littoral  drift 
which  were  computed  using  three  methods,  listed  below,  with  dredging 
data  obtained  at  the  Port  of  Visakhapatnam  on  the  east  coast  of  India. 

a.  The  Coastal  Engineering  Research  Center's  (CERC)  formula. 

b.  The  Putnam-Munk-Traylor  formula. 

c_.  Eagleson’ s formula. 

While  the  Putnam-Munk-Traylor  and  Eagleson  formulas  boch  gave  results 
that  were  within  20  percent  of  the  measured  drift,  the  amount  computed 
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using  the  CERC  formula  was  in  excess  by  a factor  greater  than  four. 
However,  all  these  methods  suffer  the  same  form  of  deficiency  as  those 
discussed  by  Galvin,  i.e.,  the  integrated  gross  transport  is  computed. 

17.  The  requirement  for  determination  of  the  velocity  field  as  a 

function  of  distance  from  the  shore  (and  therefore  a function  of  the 

water  depth)  leaves  little  choice  in  the  basic  manner  by  which  the 

longshore  velocity  profile  is  computed.  Most  studies  concerned  with 
33  3U 

this  topic  have  used  what  shall  be  referred  to  hereinafter  as 

the  radiation  stress  approach,  first  discussed  by  Lcnguet-Higgins  and 
Stewart^  and  Bowen. ^ 4 Development  of  the  radiation  stress  approach 
will  not  be  recounted  here  since  it  is  discussed  quite  eloquently  in  the 
above-referenced  papers.  A detailed  discussion  of  the  application  of 
this  approach  is  contained  in  Part  III. 

18.  Komar  and  Inman  0 compared  two  models  for  determining  litto- 
ral sand  transport  with  field  measurements  obtained  at  Silver  Strand 
Beach,  California,  and  El  Moreno  Beach,  Baja,  California.  The  first 
model  is  described  as  principally  intuitive  and  attempts  to  correlate 
the  sand  transport  rate  with  the  longshore  component  of  energy  flux 

per  unit  length  of  beach.  Komar  has  shown,  however,  that  Equation  1 can 

be  developed  under  the  basic  assumption  that  the  longshore  rate  of  trans 

port  of  sand  is  equal  to  the  longshore  component  of  sediment  carried 

forward  under  the  bore  of  the  breaking  wave.  The  second  model  was  pro- 
37 

posed  by  Bagnold  and  is  expressed  as 


U ‘ K’(ECn).  %(Vi>U» 


where 


Ie  * rate  of  transport  of  immersed  weight  of  sand  past  section  of 
beach* 

K’  * dimensionless  factor  of  proportionality 
E = wave  energy  density 
C = wave  group  velocity 


* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 
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longshore  velocity  of  unspecified  origin 


4 


I 


I 


Um  = maximum  orbital  velocity  of  waves 

= angle  of  wave  crest  at  breaking  with  respect  to  the  shoreline. 
Komar  and  Inman  concluded  that  the  models  are  equivalent  if  driving 
forces  of  the  longshore  current  other  than  those  due  to  wave  action  from 
an  oblique  direction  are  ignored. 


3.1 


Wm-Wm 
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PART  III:  OUTLINE  AND  DISCUSSION  OF  METHODOLOGY 


Introduction  and  Background 

19.  The  basic  approach  used  in  this  investigation  to  obtain  a 
quantitative  estimate  of  bed-load  sediment  transport  in  the  offshore 
zone  vas  developed  over  a period  of  about  25  years  by  the  late 

Prof.  H.  A.  Einstein  and  his  students  at  the  University  of  California, 

Berkeley.  Most  of  this  work  has  been  summarized  by  Prof.  Einstein. 

Einstein's  method  evolved  as  an  application  to  the  vave-induced  sediment 

38 

transport  problem  of  the  Einstein  bed-load  function,  which  had  origi- 
nally been  developed  to  calculate  sediment  transport  in  streams.  Al- 
though bed-load  movement  due  to  flow  in  streams  is  unidirectional, 
Einstein  and  his  students  thought  that  his  theoretical  relationships 
could  describe  the  sediment  transport  caused  by  waves.  All  of  their 
work  in  developing  this  apprcach  was  conducted  as  a series  of  laboratory 
studies  which  substantiated  the  theoretical  predictions.  The  first  step 

was  the  description  of  the  velocity  distribution  at  the  bed  due  to 
39-1*2 

waves.  The  mechanics  of  bed-load  sediment  motion  due  to  waves 

1*3  1*1* 

was  then  described  by  Kalkanis  and  later  modified  by  Abou-Seida. 

Work  on  the  mechanics  of  suspended  sediment  motion  also  has  been  under- 
taken; however,  the  precise  physical  description  of  this  process  is 
still  an  open  question  and  this  work  was  not  used  in  the  present 
study.45*16 

20.  The  model  used  for  this  investigation  will  follow  the  methods 

1*1* 

of  calculating  bed-load  transport  described  by  Abou-Seida.  It  is  as- 
sumed that  sediment  motion  is  initiated  in  the  offshore  zone  by  the 
velocity  produced  at  the  bottom  by  waves.  This  sediment  motion  is 
entirely  oscillatory  and  causes  no  net  movement  since  small-amplitude 
wave  theory  is  used.  However,  since  sediment  is  in  motion,  it  can  be 
acted  upon  by  secondary  currents  that  produce  a net  transport.  These 
secondary  currents,  such  as  longshore  and  mass  transport  currents,  are 
not  normally  intense  enough  in  the  offshore  zone  to  initiate  sediment 
motion;  however,  they  can  act  upon  sediment  already  in  motion.  In  the 
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present  study,  tnese  secondary  currents  are  described  by  the  Longuet- 

12 

Higgins  mass  transport  equation,  as  originally  recommended  by 

1*3 

Kalkanis,  and  by  the  Longuet-Higgir.s  and  Stewart  radiation-stress  anal- 
ysis of  longshore  currents. 1:> 


Summary  of  Assumptions 


21.  A variety  of  assumptions  and  limitations  were  necessary  to 

develop  a practical  method  for  the  calculation  of  bed-load  concentration: 

1*3  1*1* 

a.  The  model  developed  by  Kalkanis  and  Abou-Seida  as- 
sumes a horizontal  bed.  This  assumption  cannot  be 
precisely  satisfied  in  the  offshore  zone;  however,  bottom 
slopes  are  normally  so  gradual  that  the  gravitational 
forces  exerted  are  extremely  small  compared  with  those 
exerted  by  wave  action. 

b.  The  sediment  being  moved  can  be  described  by  its  mean 
particle  diameter.  In  most  areas  this  is  a valid  assump- 
tion since  sediment  sizes  do  not  vary  greatly  in  the  off- 
shore zone,  and  the  standard  deviation  of  most  beach 
sediment  samples  is  small. 

c_.  Sediment  particles  are  generally  spherical  in  shape. 

This  assumption  is  satisfactory  since  beach  sediment  is 
primarily  quartz  grains  that  tend  toward  a spherical 
shape  after  long  exposure  to  abrasive  action.  Also,  it 
has  been  shown  experimentally  that  the  shape  of  sediment 
particles  is  much  less  important  than  particle  size  under 
most  conditions  of  sediment  transport  .**  7 

d.  The  model  is  only  applicable  co  clastic  sediment.  No 
attempt  was  made  to  include  the  interparticle  forces 
characteristic  of  silt-  or  clay-sized  particles.  The 
model  may  be  of  questionable  accuracy  when  applied  to 
material  in  the  very  fine  sand  size  range  (3  to  It  ♦or 
0.125-  to  0.063-mm  diam).  Fortunately,  most  sediment 
found  in  the  zone  of  interest  is  normally  sand  in  the 
fine-to-mediun  grain  size  ranges. 

c_.  The  model  as  presently  designed  considers  only  bed-load 
sediment  movement.  Although  some  sediment  is  thrown 
into  suspension  in  the  offshore  zone,  it  is  not  consid- 
ered to  be  a large  enough  quantity  or  the  method  of 
describing  it  accurate  enough  to  justify  an  attempt  to 
include  it. 

£.  Einstein's  assumption  is  used  (from  his  work  with  river 
sediment  transport38)  that  the  bed  layer  is  only  two  sed- 
iment grain  diameters  thick.  Although  this  assumption 


13 


has  been  criticized  by  Brenninkmeyer  as  too  restrictive, 11 
it  is  believed  that  Einstein's  assumption  is  still  the 
best  available  for  practical  application  at  the  present 
time. 

Sfc'iment  movement  occurs  only  when  conditions  in  the 
boundary  layer  cun  be  described  as  turbulent.  Although 
this  assumption  is  questioned  by  Komar  and  Miller, ^9  who 
proposed  that  movement  can  occur  under  laminar  flow  con- 
ditions, the  turbulent  flow  criterion  establishes  oy 
Kalkanis**3  is  used. 


Outline  of  Method  Development 


22.  As  previously  noted,  the  model  outlined  herein  ts  essentially 

U3  1*1* 

that  developed  by  Kalkanis  and  modified  by  Abou-Seida.  This  brief 
introduction  is  only  a synopsis  of  the  approach  used  and  the  conclusions 
reached.  The  original  work  by  Kalkanis  or  Abou-Seida  should  be  con- 
sulted if  a more  thorough  description  of  the  methodology  is  desired. 

For  a complete  description  of  the  development  work  that  led  to  this 

it 

model,  the  summary  by  Einstein  is  r-_  commended. 

a.  Description  of  bed  motion  due  to  waves.  From  snail- 

amplitude  wave  theory  (see  Viegel^O  &nd  many  others),  the 
horizontal  particle  velocity  due  to  oscillatory  waves  at 
the  ocean  bottom  immediately  outside  the  boundary  layer 
can  be  described  by 


T sinh  (kh) 


sin  (kx  - ut) 


where 


v * velocity 

H * wave  height 

T = wave  period 

k * wave  number  * 2*/L 

h * depth  below  still-water  level 

x * distance  in  the  direction  of  wave  propagation 

<i>  - wave  frequency  * 2n/T 


By  rearranging  terms,  this  can  also  be  expressed  as 


,*  "’S'"  /'.a 


V c 


r h ii 

(vr  / 

|2  sinh  (kh)Jj 

sin  (kx  - ut) 


* (ua)  sin  (kx  - ut) 


* v sin  (kx  - ut) 
o 


(3) 


where 


a = maximum  amplitude  of  horizontal  water  displace- 
ment at  the  bottom 


v * maximum  horizontal  bottom  velocity  outside  the 
boundary  layer 

These  two  terms  (a  and  vc)  provide  a reference  from 
which  the  velocity  distribution  within  the  boundary  layer 
can  be  calculated.  The  distribution  of  velocity  within 
the  boundary  layer  is  of  the  utmost  importance  because 
all  sand  size  particles  at  rest  are  contained  in  this 
thin  layer  next  to  the  bottom. 


b.  Velocity  distributions  within  the  oscillating  boundary 
layer: 

(l)  The  velocity  distribution  within  the  laminar  boundary 
layer  can  be  described  by  Stokes  classic  linearized 
solution  to  the  Havier-Stokes  equations. 51  Using  the 
terminology  of  Kalkanis , this  is  expressed  as 


.[■ 


v |sin  ut  - e 
o 


-M 

sin  (ut  - 6 z) 
o 


i 


(*») 


where 


B « 1/6  = V^TIv 


v 

z 


boundary  layer  thickness  as  defined  by 
Schlichting52 

kinematic  viscosity 

distance  from  the  theoretical  bed 


This  solution  has  been  extended  to  the  case  of  the 
turbulent  boundary  layer  by  Kalkanis,  by  generalizing 
the  equation  for  the  laminar  velocity  distribution 
in co  the  form 

v « VQ|sin  ut  - f^z)  sin  Jut  - f2(z)]j  (5) 
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Combining  the  two  sine  functions  vector  tally,  the 
result  can  he  written  in  the  form 


v ■ v 


|l  - 2f1(z)  cos  f>(z) 

2 I1/2 

+ f*(z)j  sin  (wt  - 0) 


where 


6 * phase  angl**  * tau 


-1  T fi^  Sin  t2^  1 

[l  - t^z)  cos  fgfzjJ 


t'^(z)  * velocity  amplitude  function 

tAi)  * phase  shift  function 

<* 

Both  f^(z)  and  f2(z)  are  experimentally  deter- 
mined functions  that  take  a different  form  depending 
on  the  type  roughness  of  the  surface  over  which  the 
velocity  profile  is  being  produced.  These  empirical 
functions  were  derived  in  laboratory  experiments  for 
a smooth  bed  and  with  two-  and  three-dimensional 


roughness  and  ere  summarized  by  Einstein.  For  beds 
composed  of  exther  sand  or  Gravel,  these  functions 
are  expressed  by 


t±(z) 


-(l33z/a6  e) 

0.5e  ' ' 


f2(z)  « 0.5(8y)2/3  (7) 

The  roughness  diameter  (e)  is  assisted  equal  to  the 
mean  sediment  diameter  (D^q).  These  expressions 
were  obtained  from  experiments  that  covered  the  fol- 
lowing ranges  of  values 

0.0009  f 1 < e < 0.0453  ft 


0.10U  f t < a < 2.00  ft 


0.169  < « < 5.82 


Water  viscosity  (i.e.  temperature  and  salinity)  was 
held  fixed  during  the  experiments. 

(2)  As  a caution,  it  must  be  mentioned  that  not  all 
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investigators  agree  that  the  turbulent  boundary  layer 
velocity  profile  can  be  describe.,  in  this  Banner.  In 
his  review  of  the  current  state  of  tto  art  in  bound- 
ary layer  research,  Teleki'3*P  **  states  that  "there 
is  no  adequate  theory  for  the  mechanics  of  turbulent 
flow  at  this  time."  More  specifically,  Foda5^  noted 
that  KaAkanic*  empirical  description  of  the  turbulent 
velocity  distribution  was  not  verified  by  his  experi- 
mental data. 

£.  Flow  regime  transition  criteria.  %\nce  it  has  be<i.i:  as- 
sumed that  sediment  will  begin  to  >,ove  only  under  t.urbu- 
lent  '.md it ions , criteria  must  be  established  to  deter- 
mine what  flow  characteristics  constitute  transition  from 
laminar  to  turbulent  conditions.  The^e  criteria  have 
been  established  experimentally  firvi  the  work  of  Li  ,39 
Manohar,1*1  and  Kalkanis1*2  for  the  three  types  of  bed 
roughness  mentioned  previously.  For  beds  composed  of 
sand  and  gravel,  Abou-Seida1*1*  recommends  a Reynolds  num- 
ber criterion  of 

R « ^2.  , ioh  (8) 

c v 

For  all  Reynolds  number  values  greater  than  or  equal  to 
this  criterion,  the  flow  regime  is  considered  turbulent 
and  the  empirical  velocity  distribution  previously  de- 
scribed is  applicable.  The  criterion  Tot  the  transition 
froc  iaainar  to  turbulent  flow  conditions  is  the  subject 
of  a continuing  controversy.  In  his  review  of  the 
subject,  Teleki53  found  a wide  variation  in  the  Reynolds 
number  criteria  that  describe  this  transition  and  he 
thought  that  none  were  satisfactory.  In  particular, 
Sleath  criticises  this  Reynolds  number  criterion  for 
being  based  on  insufficient  data. 55 

d.  Motion  of  sediment  due  to  wave  notion: 

(1)  Once  the  velocity  distribution  within  the  turbulent 
boundary  layer  is  known,  the  amount  of  sediment  in 
motion  can  be  found  by  use  of  the  bed-load  equation 
derived  by  Kalkanis**3  and  extended  by  Abou-Seida.1*1* 
This  derivation  follows  the  probabilistic  procedure 
developed  by  Einstein. 38  since  a derivation  of  this 
equation  is  rather  lengthy,  only  most  applicable 
information  will  be  described  here.  Einstein  pro- 
vides a short  synopsis  of  the  derivation.1* 

(2)  The  bed-load  equation  establishes  a /elation  between 
two  nondimensional  parameters  that  describe  the  fluid 
system  and  the  resulting  sediment  transport.  The 
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first  of  these,  the  flow  intensity  function  is  given 
as 


where 

i * flow  intensity  function 

* density  of  the  sediment  particle 

Pf  * density  of  the  fluLc. 

g = acceleration  due  to  gravity 

D50  * mean  sediment  diameter 

The  velocity  va  is  the  mean  unsteady  velocity  within 
the  boundary  layer  given  by  the  equation  for  the  tur- 
bulent velocity  distribution  at  a point  0.35  D^q  above 
the  theoretical  bed.  The  velocity  at  this  distance 
was  determined  by  Einstein  and  El  Samni^b  as  best  de- 
scribing the  lift  forces  on  each  sediment  grain  for 
unidirectional  flows.  Whether  this  applies  equally 
as  well  for  oscillatory  notion  has  not  yet  been  deter- 
mined, but  it  is  a reasonable  assumption. 

(3)  The  second  of  these  nondimensional  equations  is  the 
bed-load  function  given  by 


where 

qn  = oscillatory  bed-load  rate 

Yg  * unit  dry  weight  of  bed  material 

(M  The  relation  between  the  flow  intensity  function  and 
the  bed-load  function  is  given  by  the  bed-load  equa- 
tion ( Abou-Seida^ ) 
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where 


|p%-  *. 


m ~ a rundom  variable  that  has  a nomal  distri- 
bution of  zero  mean  and  a standard  devia- 
tion of  one 

The  terms  A , B , find  Hq  are  empirically  derived 
constants  which  best  fit  the  theoretical  curve  to  the 
experimental  data  of  Kalkanis  and  Abou-Seida  (see 
Figure  1).  The  best  fit  was  obtained  with  values  for 
these  terms  of: 


A = 13.3 
B * 6.0 

n = jo 
o 

(5)  The  flow  intensity  function  as  expressed  in  the  bed- 
load equation  has  been  modified  by  the  following 
expression  to  account  for  a particular  physical 
phenomenon: 

= « (12) 

The  term  t is  a correction  factor  introduced  by 
Einstein^  and  commonly  known  as  the  "hiding  factor.” 
This  factor  was  introduced  to  account  for  the  fact 
that  under  turbulent  flow  conditions  very  small  sedi- 
ment particles  are  partially  or  completely  submerged 
in  the  laainnr  sublayer  where  they  are  not  affected 
by  the  velocity  fluctuations  in  the  fully  turbulent 
.lone.  For  the  case  of  oscillatory  notion,  Abou-Seida 
has  determined  this  correction  experimentally.  He 
has  found  that  the  correction  is  related  to  the  param- 
eter D5o/i.39  , where  5*  is  the  thickness  of 

the  laminar  sublayer.  From  dimensional  considera- 
tions and  experimental  results,  Einstein  and  Li^O 
found  that  the  thickness  of  the  laminar  sublayer  in 
oscillatory  flows  can  be  expressed  as 


o 


Using  this  expression,  the  parameter  D50/I.39  5* 
can  be  computed  and  the  hiding  factor  can  be  deter- 
mined from  the  plot  of  Abou-Seida' s experimental 
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Figure  1.  Plot  of  bed-load  equation;  data  from  Kalkanis  J and  Abou-Seida 


results  (Figure  2)  or  maerically  as  vill  be  illus- 
trated later. 

(6)  Values  of  the  bed-iced  function  frost  the  corrected 
flov  intensity  function  can  be  obtained  by  use  of  the 
curve  developed  from  experimental  data  (see  Figure  1) 
or  by  a numerical  solution  of  the  bed-load  equation 
that  vill  be  described  later. 


(?)  When  a value  of  the  bed-load  function  is  obtained,  a 
value  of  the  oscillatory  bed-load  rate  qB  also  can 
be  established.  This  represents  the  average  amount 
of  sediment  in  motion  due  strictly  to  the  oscillatory 
wave  action;  qB  is  expressed  in  terms  of  the  dry 
weight  of  sediment  in  motion  per  unit  of  time  per  unit 
vidth  of  beach  normal  to  the  direction  of  wave  prop- 
agation. As  such,  qs  represents  no  net  sediment 
transport. 

(8)  Kalkanis**3  obtained  a value  of  net  sediment  transport 
in  the  direction  of  the  secondary  current  by  first 
reducing  q0  to  an  average  bed-load  concentration 
by  developing  the  following  expression 


M 


B 


w 


2D_  v 
50  n 


where 


Cg  * average  bed-load  concentration 
M * coefficient  of  proportionality 
2D50  * thickness  of  the  bed  layer  (Einstein^®) 
average  velocity  at  a distance  D50 


m 


from  the  theoretical  bed 


In  this  expression,  the  average  velocity  is  calcu- 
lated from  the  equation  for  the  turbulent  velocity 
distribution  at  a distance  of  D^q  from  the  theoret- 
ical bed.  The  coefficient  M has  been  theoretically 
evaluated  by  Kalkanis  for  small-amplitude  waves  to 
have  a value  of  0.6l8.  Abou-Seida  found,  however, 
that  for  steep  waves  with  short  periods  and  wave- 
lengths, the  value  of  M must  be  larger  to  account 
for  the  fact  that  acceleration  effects  ignored  in 
Einstein's  original  development  of  the  bed-load  equa- 
tion can  no  longer  be  neglected.  At  first,  it  may 
not  seem  logical  that  a correction  for  acceleration 
effects  should  be  applied  to  the  resulting  sedine:  t 
transport  rather  than  to  the  flow  intensity  func- 
tion. However,  as  Einstein  has  pointed  out the 
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acceleration  force  applies  not  only  to  the  fixed  sedi- 
ment particles  on  the  bed  but  to  those  in  motion  as 
well. 

(9)  Using  experimental  data  of  his  own  in  addition  to 
those  of  Vincent, 57  Abou-Seida  plotted  values  of  M 
against  the  parameter  a/Djg  » which  he  concluded  was 
an  appropriate  indicator  of  the  acceleration  effects 
(Figure  3).  This  shows  that  for  values  of  a/D^g 
less  than  about  300,  the  coefficient  M increases 
sharply  as  an  ix  dication  of  increased  acceleration 
effects.  For  values  of  a/Dcg  greater  than  300  s 
values  of  M remain  remarkably  constant  at  about  the 
theoretical  value  of  0.6l8.  For  most  conditions  in 
the  offshore  zone  where  small  amplitude  theory 
applies,  the  theoretical  value  is  satisfactory  since 
acceleration  effects  may  be  neglected. 

e.  Net  bed-load  sediment  transport; 

(1)  The  average  concentration  of  bed-load  material  (Cg) 
due  to  oscillatory  wave  action  provides  a measure  of 
the  amount  of  sediment  in  motion.  This  is  not  a mea- 
sure of  sediment  transport  since  it  is  due  to  the 
periodic  wave  motion  that  develops  no  net  movement 
(within  the  approximations  of  the  Airy  wave  theory). 
However,  since  this  material  is  already  in  motion,  it 
can  be  moved  readily  by  any  translatory  current. 

Such  currents  can  be  tidal,  or,  as  is  more  prevalent 
in  the  offshore  zone,  currents  developed  by  wave 
action  such  as  the  mass  transport  and  longshore  cur- 
rents. By  superposing  such  a current  on  the  average 
bed-load  concentration,  the  resulting  net  sediment 
transport  can  be  expressed  by  the  following  equation 


(15) 


where 


= net  bed-load  sediment  transport 


U(z)  = depth-dependent  velocity  function 


This  equation  gives  a rediment  transport  in  the 
direction  of  the  velocity  vector.  Qg  is  expressed 
as  the  dry  weight  of  sediment  transported  per  unit 
time  per  unit  length  normal  to  the  direction  of  the 
velocity  vector. 

(2)  Kalkanis  originally  recommended  that  the  mass  trans- 
port current  as  it  is  described  by  Longuet-Higgins-*-^ 
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be  applied  as  the  tranalatory  current  to  produce  net 
sediment  transport.  As  a matter  of  interest,  it  was 
from  this  application  of  the  mass  transport  velocity 
to  his  laboratory  experiments  that  Abou-Seida  de- 
veloped most  of  the  correction  factors  to  the  bed- 
load equation.  By  simply  modifying  the  notation  of 
the  Longuet-Higgins  equation,  Kalkanis  expressed  the 
mass  transport  velocity  as 

„2  . / -B  z -2B  z\ 

U(z)  * r 15  - 8e  0 cos  Bz  + 3e  0 ) (l6) 

16  sinn  kh  \ / 

Net  sediment  transport  in  the  direction  of  wave  prop- 
agation can  be  calculated  by  using  this  expression  in 
conjunction  with  the  equation  for  net  transport. 

Description  of  Longshore  Velocity  Field 

23.  The  longshore  velocity  field  is  calculated  from  that  described 

5ft 

by  Longuet-Higgins.  The  assumptions  inherent  in  this  calculation  are: 

a.  The  longshore  velocity  field  is  stationary. 

b.  The  beach  is  planar  and  of  infinite  length. 

c_.  Velocity  is  independent  of  depth  in  the  y direction  but 
return  flow  is  allowed  in  the  x direction. 

d.  Currents  due  to  other  than  wave  action  (e.g.,  wind- 
driven  currents)  are  sufficiently  small  that  their  inter- 
action with  the  wave  regime  is  inconsequential. 

e.  Coriolis  force  5°  neglected. 

f_.  Sediment  size  is  uniform  over  the  beach. 

g.  Beach  slope  is  sufficiently  small  so  that  the  effects  of 
wave  reflection  are  negligible. 

Assumption  (3)  does  not  permit  any  return  flow  in  the  form  of  rip  cur- 
rents. Evidence  indicates  that  the  return  of  the  mass  flow  due  to  wave 

advance  takes  place  at  intermediate  depths.  This  has  been  discussed  in 
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detail  on  a theoretical  basis  by  Longuet-Higgins. 

24.  No  attempt  was  made  during  this  study  to  directly  verify  the 

validity  of  the  radiation  stress  approach.  Bowen  conducted  extensive 

laboratory  experiments  and  found  good  correlation  between  the  predicted 

and  observed  values  of  the  water-surface  setdovn  over  the  range  of  wave 
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heights  and  beach  slopes  tested.  This  agreement  presumabJy  infers 
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that  confidence  can  be  placed  in  the  velocity  field  obtained  by  this 
method. 

25.  The  basic  equations  used  in  this  development  are  those  which 
appear  in  References  1 6 and  58.  The  equations  of  motion  used  in  the 
radiation  stress  approach  are  (see  Figure  k for  definition  of  coordinate 
system) 


h (y  f)  - ■ 0 


where 


t * net  stress  in  y direction 
y 

Ue  = coefficient  of  horizontal  eddy  viscosity 
3 \ » mean  stress  on  bottom 


a mean  stn 
and,  in  shallow  water 


i e2Pf(B»)3/2s(y-2)  , x < x. 


(By)  » f 6Spf(gh)1/2v 


where  S is  a friction  coefficient  (see  paragraph  26). 
eddy  viscosity  is  calculated  from 


The  horizontal 


ne  = Npfx(gh)J 


where  N is  a coefficient  of  proportionality.  From  the  expected  maxi- 
mum turbulent  velocities  and  maximum  mixing  length  Longuet-Higgins  rea- 
sons that 

0 < N < 0.016' 


Assuming  a constant  beach  slope 
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where  s * dh/cx  . Equation  17  can  then  be  rewritten  as 

l 3/2 

JU  / 5/2  a»\  1/2  , }-r*  • 

’•*(*  •*)"*  jo 

where  the  values  of  p , q , and  r are  given  as 


P - «V1/2s3/2 
1 * f 6SPi*1/2«1/2 

..Si  3/2  5/2  ,ln  “b 

' f6P**  * ^ 


The  nondimensional  variables 


y = i_  v = 2L_ 
% * *. 


are  introduced  where  v is  defined  as 

o 


51  i ^/2, 


o 8 S 


(ghb)J 


"s  sin 


Frosi  Equations  21  and  22 


2 fe  (*5/2  i)  - ^ 


-x3/2  , 0 < X < 1 
0 , 1 < X < «• 


where 
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\2Aes/ 


P is  a measure  of  the  horizontal  nixing.  The  boundary  conditions 
imposed  upon  . at  the  breaker  line  (X  * l)  are  that  both  V and 


3V/3X  are  continuous.  The  complete  solution  to  Equation  2 f is 


X + AX  , 0<X<1 


, 1 < X < - 


(28) 


vhere  the  values  of  A , , and  are  given  by 


P d — 

’ ^5 


px  ■ - f + T p) 
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(29) 


(30) 


p2  - 1 

B * — A 

1 pl-p2 


B„  * 


Pi-1 


2 Pl-P2 


(31) 


The  nondimens ional  form  of  the  velocity  profile  is  obtained  by  substi- 
tuting Equations  27*  29,  30,  and  31  into  Equation  28.  Equations  2k 
and  25  are  used  to  relate  results  of  Equation  28  to  the  prototype 
problem.  Figure  5 shows  a family  of  current  profiles  as  a function  of 
the  mixing  parameter  P . This  figure  illustrates  the  relative  insensi- 
tivity of  the  velocity  profiles  to  the  value  of  P (V  varies  less  than 
an  order  of  magnitude  with  values  of  P over  five  orders  of  magnitude). 

26.  The  calculation  of  typical  values  of  the  necessary  parameters 
is  illustrated  as  follows.  Several  preliminary  computations  were  made 
using  values  of  N such  that  0.005  < N < 0.15  . Results  of  these 
calculations  were  examined,  and  it  was  concluded  that  the  system  is  rel- 
atively insensitive  to  the  values  of  N within  the  indicated  range. 

N was  assigned  the  value  of  0.01  for  the  computations  made  for  this 
study.  An  order  of  magnitude  estimate  of  the  Reynolds  number  was  needed 
to  determine  a value  of  the  frictional  coefficient  S . Using  the 

eg 

values  in  Equation  k5  of  Longuet-Higgins? 
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Figure  5.  The  form  of  longshore  current  profiles  as  given  by 
Equation  28  sequence  of  values  of  the  mixing  parameter  P . 
Adapted  from  Longuet-Higgins5i 


where 


$ » 0.1* 

r 

g * 10  m/sec* 


u * 0.8  radian/ sec 
v « 1.3  x 10  m^/sec 


gives  the  result 


Jifgh 


R * 3 X 106  - 0(6) 


The  above  values  correspond  to  a wave  with  a 5-sec  period  and  2-m  break- 
ing height.  The  second  parameter  needed  is  the  ratio  a/e  where  e is 
a typical  value  for  a roughness  element  (for  our  purposes  c **  0.25  cm) 
and  a is  the  horizontal  excursion  of  a water  particle  given  by 


thus,  for  this  example. 


a 


>10  x 2 

~T” 


5.6  m 


and 
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0.0025 


2.2  x io 
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Reading  from  Figure  6 (from  Prandtl^ ) , 


these  values  correspond  to  a 


<*/£  = 

10* 


10* 


10* 


10* 


Figure  6.  Values  of  the  friction  coefficient  S for  flow 
over  rough  plates  (from  Prandtl6l) 

friction  coefficient  S * 0.008  . If  the  parameters  N , S , any  £ 
are  fixed,  P is  a function  of  the  bottom  slope. 

27.  To  refract  and  shoal  the  deepwater  waves,  the  governing  equa 
tions  are  those  typically  used: 


sin  a 


(35 


Snell's  law 


Dispersion 

relation 


— * tanh  Eh 


136) 


Conservation  of 
energy  flux 


H_ 

H 


*r“s“fp 


(37) 


vhere  H , 5 , and  a are  functions  of  x . From  Thornton,1**  if  per- 

colation is  ignored,  the  energy  dissipation  due  to  friction  is  given  by 


J£_  / uM  \3 

nf  ~ 12»  ysinh  Kd  / 


(38) 


vhere  f is  obtained  from 
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JL— 
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V U J u 


-0.0S  + an  — 
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and 


(39) 


v * H i — - 

Td  2 sinh  Kd 


(UO) 


Equation  39  is  solved  by  an  iterative  procedure. 

28.  To  summarize  the  procedure  thus  far:  the  deepvater  wave 

regime  is  refracted  and  shoaled  into  shallow  water  to  the  point  of  break- 
ing, taking  account  of  the  effect  of  dissipation  due  to  bottom  friction. 
Distance  from  shore,  angle  of  breaking,  and  breaking  wave  height  are 
used  to  compute  the  longshore  current  profile,  following  the  procedure 
of  Longuet -Higgins  previously  outlined.  During  the  numerical  shoaling 
and  refraction  procedure,  the  values  of  water  depth,  wave  height,  wave- 
length, and  refraction  angle  are  retained  for  use  in  computation  of  the 
bed  load. 
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PART  IV:  DffLEMBTEATIOH 

29.  To  obtain  a measure  of  the  distance  rate  of  decrease  of 
velocity  within  the  boundary  layer,  it  ***  assumed  the  boundary  layer 

co 

velocity  profile  is  represented  by  the  empirical  formula  (Schlichting^ ) 


where  Ua  is  the  fluid  velocity  outside  the  boundary  layer  and  n is 
an  integer  determined  by  the  magnitude  of  the  Reynolds  number.  For 
R ~ 10^  (turbulent). 


u ' Mi 


where  6 is  thickness  of  the  boundary  layer.  Because  of  the  difficulty 
in  determining  this  distance,  it  is  usually  taken  as  the  distance  above 
the  solid  boundary  where  the  velocity  reaches  99  percent  of  the  velocity 
of  the  well-developed  flow  exterior  to  the  boundary  layer.  This  latter 
thickness  is  denoted  by  6*  (from  Abou-Seida  ) , 


llOOv 

2U_ 


and  the  boundary  condition 


Ufi*  * j-  at  z * 5* 


Combining  Equations  1*2  and  U3  and  applying  the  boundary  condition 


/ 1100v\ 
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thence  solving  for  5 , 


2 x llOOv 
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—5  2 

For  a kinematic  viscosity  of  1 * 10  ft  /sec 


0.352 


Figure  7 is  a plot  of  the  boundary  layer  thickness  as  a function  of  the 
far-field  fluid  velocity.  For  the  purposes  of  our  discussion,  U is 

0D 

the  longshore  current  velocity  obtained  using  the  radiation-stress  ap- 
proach. The  velocity  within  the  boundary  layer  is  designated  by 
where 


UT  ■ u-(!)  . ‘-Vs 


i 


111 


It  was  decided  to  try  the  integrated  mean  value  of  U , , in  a 

first  attempt  to  compute  the  velocity  within  the  turbulent  boundary 
52 

layer.  Schlichting  gives 


UT  = 0.79  Ua 


for  n * 6 . Justification  for  this  can  be  found  by  examining  Figure  7 
where  it  is  evident  that  U reaches  a value  of  about  0.7U  within 

i ** 

the  first  1/10  of  the  boundary  layer  thickness. 

30.  Three  types  of  information  are  needed  to  calculate  an  average 
bed-load  concentration.  These  are  the  characteristic  wave  parameters 
(period,  wavelength,  height,  and  depth  of  water),  the  fluid  character- 

5 sties  (density  and  viscosity),  and  the  sediment  characteristics  (den- 
sity, mean  particle  diameter,  and  dry  specific  weight).  The  sediment 
and  water  characteristics  are  known  initially  and  wave  parameters  can  V 
calculated. 

31.  After  the  values  of  a few  necessary  parameters  are  obtained 

(a  and  u ),  the  first  step  in  the  computation  is  to  determine  if  the 
o 

requirement  of  turbulent  flow  conditions  at  the  bed  is  satisfied.  This 
is  simply  accomplished  by  checking  to  see  if  the  Reynolds  number  cri- 
terion for  turbulent  conditions  is  satisfied 
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If  this  requirement  is  met,  the  next  step  is  to  calculate  the  average 
velocity  within  the  boundary  layer  at  a distance  of  0.35  from  the 

theoretical  bed.  This  is  accomplished  by  using  the  following  modified 
form  of  the  turbulent  velocity  distribution  equation  and  the  recommended, 
forms  of  the  velocity  amplitude  and  phase  shift  functions. 

1 1/2, 


,[l  - 21'^a)  cos  f2(a)  + f2(a)]  (f) 


f-jU)  * 0.5e 


-(fc6.55/aBo) 


f2(a)  = 0.2U8  (BoD50): 


2/3 


(50) 


From  this  result,  the  flow  intensity  function,  <l>  , is  calculated  using 
Equation  9.  The  hiding  factor  l,  must  be  determined  from  a statistical 
evaluation  of  Abou-Seida's  original  data  (Figure  2).  Upon  close  exam- 
ination of  these  data,  four  distinct  groups  of  points  were  observed.  A 
least-squares  analysis  was  performed  for  each  group.  This  provided  the 
following  empirical  estimate  for  the  hiding  factor: 
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1.396 
for 


X > 0.14230 


C « 0.7315X~°*363U  for  0.4230  > X > 0.2173 


C « 0.3701X"0,8097  for  0.2173  > X > 0.1379 


5 « 0.01407X”1,921*1 


for 


X < 0.1379 


32.  Once  the  hiding  factor  is  found,  the  modified  flow  intensity 
function  is  computed  using  Equation  12  and  can  then  he  used  to  solve  the 
bed-load  equation.  Abou-Seida  integrated  the  right-hand  side  of  the 
bed-load  equation  by  using  the  tables  of  the  normal  probability  integral 
in  conjunction  with 
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To  evaluate  this  set  of  equations,  Abou-Seida  set  N equal  to  9 so  that 
represents  the  midpoint  of  one  of  the  nine  equal  intervals  dividing 

X 

the  quarter  of  the  cycle.  For  the  present  ipplication,  a power  series 
solution  is  used  to  compute  values  of  normal  probabilities.  A value  of 
90  is  used  for  H in  calculating  the  bed-load  function. 

33.  The  oscillatory  bed  load  is  given  by 


% m V 


(»> 


and  the  bed-load  concentration  is  calculated  from  Equation  lh.  v is 

x& 

found  by  calculating  the  average  velocity  within  the  turbulent  boundary 
layer  at  a distance  above  the  theoretical  bed.  M is  determined 

from  a least-squares  evaluation  of  Vincent’s^  and  Abou-Seida* s^  orig- 
inal data  in  conjunction  with  the  theoretical  value  derived  by 


f^spupsiB**' 
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Kalkanis.  This  provides  the  following  formulas  for  M 


Z = 


a 


50 


(55) 


M * 0.6180 

for  Z > 256.6333 

M = 251.905SZ-1*0831* 

for  Z <.256.6333 

With  a value  of  the  bed-load  concentration  Cg  , net  transport  can  now 
be  calculated  for  any  secondary  current  using  Equation  15.  For  the 
mass  transport-  velocity,  evaluation  of  Equation  15  can  then  be  accom- 
plished by  a simple  integration  of  Equation  16. 
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PART  V:  MODEL  VERIFICATION 


3k.  Verification  of  any  model  of  coastal  processes  is  extremely 
difficult  due  to  the  scarcity  of  reliable  and  quantitative  prototype 
data.  In  this  particular  case,  information  on  the  onshore  transport  of 
sediment  is  practically  nonexistent.  However,  estimates  of  the  net  long- 
shore transport  of  sediment  are  numerous  and  are  based  on  long-term 
studies  that  increase  their  averaged  accuracy.  It  was  therefore  decided 
to  verify  the  present  model  by  using  it  to  calculate  longshore  transport 
for  a section  of  coast  that  has  a relatively  well-established  sediment 
budget  together  with  reasonable  estimates  of  longshore  sediment  trans- 
port. Point  Pederneles,  California,  a rocky  promontory  located  approxi- 
mately 2 miles  north  of  Point  Arguello,  was  selected  for  the  following 


reasons : 

a.  The  immediate  offshore  bathymetry,  characteristics  of  the 
bottom  sediment,  mechanisms  of  sediment  transport  in  the 
offshore  zone,  and  sediment  budget  for  this  section  of 
coast  have  been  well-documented. 8,60,62 

b_.  A deepwater  wave  hindcast  station  has  been  established 
immediately  offshore  from  this  location  (31»°30,N, 

121°00’W)  and  is  based  on  data  for  three  consecutive 
years — 195^ » 1957,  and  1958  (National  Marine 

Consultants^ ) . 

c_.  Point  Pedernales  is  a rocky  headland  that  protrudes  well 
into  the  surf  zone,  thus  effectively  blocking  any  long- 
shore '.ransport  of  sediment  in  the  surf  zone.  This 
creates  a situation  in  which  the  limitation  of  the  model 
to  the  offshore  zone  does  not  prohibit  its  use  in  deter- 
mination of  the  total  longshore  transport. 

Sources  and  Quality  of  Input  Data  Used 

Wave  data 

35.  The  wave  data  used  for  the  computation  of  sediment  transport 
at  Point  Pedernales  are  those  presented  for  Station  5,  located  at  3**.5°N, 

go 

121. 0°W,  of  the  National  Marine  Consultants.  The  significant  wave 

heights  tabulated  in  the  hindcast  were  converted  to  root-mean-square 
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heights  by  dividing  by  the  factor  1 . Ul6  (cf.  Longuet-Higgins  ) . It 
is  thought  that  the  root-mean-square  wave  height  was  a more  accurate 
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measure  of  the  energy  of  the  wave  regime  than  the  significant  wave 
height,  and  this  accuracy  would  be  reflected  in  the  computation  of  the 
longshore  transport.  Since  the  coastline  in  this  area  runs  in  approxi- 
mately a straight  line  from  north-northeast  to  south-southwest,  only 
waves  from  directions  of  north  to  southwest  are  converging  on  the  study 
area  (Figure  8).  Waves  (combined  sea  suad  swell)  from  north  and  south- 
west higher  than  h ft,  the  highest  listed  as  12  ft,  occur  a combined 
total  of  less  than  one  percent  of  the  year  and  sure  therefore  judged  to 
make  an  insignificant  contribution  to  the  net  annual  transport.  Waves 
from  west-northwest  are  approximately  normal  to  the  coast  and  although 
they  would  produce  an  onshore  transport  of  sediment,  they  would  not  con- 
tribute to  the  total  longshore  transport.  Therefore,  only  waves  from 
north-northwest,  northwest,  west,  and  west-southwest  were  used  in  the 
longshore  sediment  transport  computations. 

36.  By  averaging  wave  directions  into  arcs  of  22.5°,  the  hindcast 
introduces  the  single  most  important  source  of  error  into  the  model. 

The  longshore  velocity  and  it»  associated  sediment  transport  sure  very 
sensitive  to  the  angle  of  wave  incidence  and  the  amount  of  error  thus 
introduced  can  be  significant. 

Fluid  data 

37*  Basic  physical  properties  of  seawater  were  obtained  for 
nearby  Avila  Beach,  California  (National  Ocean  Survey^),  from  daily 
observations  made  at  the  tide  recording  station.  The  following  annual 
averages  were  used  in  the  computations. 

Temperature  (T)  * 56.8°F  = 13.8°C 
Density  (pu)  = 1.0252 
Salinity  (S)  * 33.6  °/oo 

The  value  of  salinity  was  converted  to  chlorinity  by  use  of  the  empir- 
ical relationship  recommended  by  Sverdrup,  Johnson,  and  Fleming.^’*’  ^ 
Although  more  precise  formulas  have  been  developed,  this  simple  relation 


Uo 


is  stiff icient  for  the  present  application.  The  resulting  chlorinity 

vast 


Chlorinity  ( CL)  = 18.6  °/oo 

From  these  data  and  the  relationship  developed  by  Miyake  and  Koizumi,^* 
the  dynamic  viscosity  is: 

U = 1.26(10  ) poise 

Using  the  established  average  density,  the  kinematic  viscosity  is: 

v « 1.23(10“2)  stokes  = 1.32(10-5)  ft2/sec 

The  value  of  viscosity  used  can  be  a source  of  minor  error  since  vis- 
cosity is  sensitive  to  temperature  and  also  varies  with  salinity.  Er- 
rors also  can  be  introduced  In  areas  with  spatial  variations  in  tempera- 
ture, density,  and  salinity.  Only  slight  variations  of  these  parameters 
are  present  in  the  area  under  investigation. 

38.  Preliminary  calculations  with  the  model  indicate  that  sedi- 
ment transport  increases  with  water  temperature.  This  is  similar  to  the 
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trends  found  by  Franco0  and  by  Taylor  and  Vanoni00  for  unidirectional 
flows.  The  model  also  indicates  a substantial  increase  in  Che  depth  at 
which  sediment  movement  is  initiated  with  increasing  temperature.  This 
is  logical  since  the  increasing  water  temperature  lowers  the  viscosity, 
thus  reducing  the  magnitude  of  the  velocity  necessary  to  induce  turbu- 
lent flow. 

Bathymetric  data 

39.  One  of  the  assumptions  presently  contained  in  the  model  is 
that  the  bathymetry  is  represented  by  a simple  plane  beach.  Therefore, 
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the  data  given  by  Trask  and  Bower  and  Inman,  and  National  Ocean  Sur- 
vey Hydrographic  Surveys  H-57^6,  H-57^1,  and  H-5799  were  used  to  develop 
an  average  offshore  slope  of  0.0181*.  Figure  9 illustrates  the  slope 
used  in  the  model  as  compared  with  the  actual  slope  off  Point  Pedernoles 


1*2 


in  the  immediate  offshore  zone.  The  use  of  an  averaged  slope  in  this 
area  can  be  expected  to  result  in  a higher  conputed  transport  rate  than 

g 

the  use  of  the  actual  slope  because  Trask0  believes  all  sediment  move- 
ment in  this  area  takes  place  in  depths  shallower  than  60  ft. 

Sediment  data 

1*0.  The  assumption  that  the  model  predicts  the  sediment  transport 
on  a beach  of  a constant  uniform  grain  size  coincides  well  with  actual 

O 

field  data  available  for  Point  Pedernales.  Trask  shows  that  for  the 
entire  offshore  area  shown  in  Figure  9»  the  mean  grain  size  of  all  sedi- 
ment samples  taken  varied  from  only  1 h6  to  153  microns;  therefore,  an 
average  grain  size  of  0.15  mm  (or  U.92  * 10-1*  ft)  was  used. 

1*1.  It  was  assumed  that  the  sediment  was  predominantly  quartz 

sand  with  a density  of  2.65.  It  was  also  assumed  that  the  dry  unit 
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weight  of  the  sediment  was  109  pcf  as  given  by  Terzaghi  and  Peck  for 
dense,  uniform  sand  although  values  as  low  as  90  pcf  can  be  expected  for 
loose,  ’ini form  sand.  The  larger  value  was  used  because  saturated  fine- 
grained sand  tends  to  be  densely  packed. 

Results  of  Verification  Computations 

•*2.  Using  the  model  as  presently  developed,  annual  longshore 
sediment  transport  volumes  were  calculated  for  every  applicable  combina- 
tion of  intervals  of  wave  period,  wave  height,  and  incident  angle  with 
all  other  conditions  remaining  constant.  These  values  were  than  totaled 
for  both  upcoast  arid  downcoast  sediment  transport,  and  the  gross  and  net 
transport  was  computed.  Two  separate  computations  were  conducted  in 
this  manner.  The  first  computation  considered  all  sediment  transported 
in  the  offshore  zone  (i.e.,  outside  the  breaker  zone),  while  the  second 
considered  only  transport  in  the  offshore  zone  in  depths  greater  than 
12  ft.  This  latter  computation  was  conducted  due  to  the  precipitous 
manner  in  which  Point  fcdernales  drops  into  the  Pacific  Ocean,  and  its 
general  rocky  nature.  Analysis  of  the  bathymetric  data  shows  only 
depths  of  18  and  12  ft  in  close  proximity  to  the  shore.  The  depth  of 
12  ft  was  arbitrarily  established  as  the  base  of  the  rocky  promontory 
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and  the  beginning  of  the  sediment-covered  slope.  Results  of  these  com- 
putations are  shown  below  along  with  the  longshore  transport  estimates 
of  Bowen  and  Inman^0  for  Surf,  California,  located  1*  miles  north  of 
Point  Pedernales. 
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Longshore  Transport  Estimates  for  Point  Pedernales  (in  10  ydJ/yr) 


Reference 

U pc oast 
Transport 

Downcoast 

Transport 

Gross 

Transport 

Net 

Transport 

Downcoast 

Bowen  and  Inman**® 

251* 

307 

561 

53 

Model  (all  offshore 
zone  transport) 

305 

519 

821* 

2ll* 

Model  (offshore  zone 
transport  in  depths 
greater  than  12  ft) 

217 

2l*0 
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23 

Discussion  of  the  Results  of  Verification 

U3.  The  above  tabulation  shows  that  an  approximate  verification 
was  achieved  when  the  results  of  the  computations  were  considered  in 
relation  to  the  assumptions  made.  Reliability  of  the  data  used  for 
verification  is  somewhat  unknown  although  it  is  certainly  as  good  as  or 
better  then  any  other  such  data.  The  longshore  transport  estimates  of 
Bowen  and  Inman*’®  are  based  on  calculations  of  the  longshore  component 
of  wave  power.  Although  this  method  seems  to  work  out  quite  well  in  the 
overall  sediment  budget  for  this  particular  region  of  the  California 
coast,  it  is  not  based  on  observed  prototype  measurements;  it  is,  how- 
ever, the  best  information  available. 

1*1*.  Computations  made  with  the  model,  considering  all  longshore 
transport  in  the  offshore  zone,  show  consistently  high  results  relative 
to  the  estimates  of  Bowen  and  Inman.  This  could  be  because  no  shallow- 
water  cutoff  point  for  sediment  transport  was  established  for  this  com- 
putation or  it  could  be  due  to  some  of  the  assumptions  inherent  in  the 
method  of  calculation.  It  should  be  reemphasized  that  the  model  assumes 
sediment  is  available  for  movement. 

1*5.  A much  better  estimate  of  the  sediment  transport  is  achieved 
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by  restricting  the  transport  to  depths  in  excess  of  12  ft;  justification 
for  this  restriction  is  based  on  the  existence  of  the  rock  promontory. 
Although  smaller  than  the  estimate  of  Bowen  and  Inman,  the  calculated 
values  for  northerly  and  southerly  transport  are  similarly  proportioned. 
Direction  of  the  net  drift  is  the  same  as  for  the  estimate  which  in- 
cludes movement  in  depths  less  than  12  ft;  however,  its  relative  magni- 
tude is  significantly  reduced.  This  is  primarily  due  to  the  fact  that 
the  waves  coming  from  the  northwest  and  north-northwest  sure  of  rela- 
tively low  heights  for  most  of  the  year,  and  will  produce  sediment 
movement  only  in  shallower  depths.  Waves  from  the  west-southwest  and 
vest  are  higher  and  produce  sediment  movement  in  deeper  water.  There- 
fore, by  establishing  a 12-ft  cutoff  point  as  the  shallowest  depth  at 
which  sediment  can  move  around  the  promontory,  waves  from  the  west- 
southwest  and  west  produce  a greater  relative  effect  on  the  estimate. 
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PART  VI:  RECOWtEKDATIONS 


U6.  Recommendations  are  of  two  general  categories:  those  of  a 

general  nature  pertaining  to  the  verification  of  the  model  and  those 
that  are  more  specific  pertaining  to  modifications  of  the  model. 

a.  With  few  exceptions,  most  sediment  transport  models  in 
the  past  have  used  previously  existing  wave  data  of  un- 
known or  questionable  reliability  or  gathered  wave  data 
on  a "make-do”  basis.  Therefore,  it  is  recommended  that 
a coordinated,  well-planned  program  be  undertaken  to 
gather  wave  data.  For  example,  a concerted  effort  must 
be  made  when  designing  and  installing  wave  data  instru- 
mentation to  ensure  that  it  is  capable  of  performing 
properly  during  periods  of  extremely  severe  wave  action, 
since  it  is  during  these  periods  that  the  potential  for 
sediment  movement  is  greatest.  Concurrently,  detailed 
bathymetric  surveys  of  the  verification  area  should  be 
performed.  In  the  event  of  a storm  severe  enough  to 
modify  the  morphology  of  the  immediate  area  in  a few  days, 
the  usual  sediment  tracing  techniques  might  pr^ve  inade- 
quate to  provide  information  on  the  movement  of  such 
great  amounts  of  sediment.  One  could  then  resort  to  a 
poststorm  survey  to  obtain  a gross  measure  of  the  magri- 
tude  of  the  movement.  Moreover,  it  is  important  that  the 
distribution  of  sediment  sizes  be  determined  in  order  to 
verify  a model.  This  could  be  done  in  conjunction  with  a 
bathymetric  survey  and  also  would  provide  valuable  infor- 
mation on  the  redistribution  of  sediment  sizes  in  the 
event  of  a severe  storm.  An  array  of  current  meters 
should  be  installed.  In  addition  to  meters  positioned  to 
determine  current  velocity  as  a function  of  distance  from 
shore,  there  should  be  at  least  one  location  near  the 
breaker  zone  where  the  vertical  velocity  distribution  is 
measured. 

b.  A list  of  recommended  modifications  in  order  of  their 
estimated  overall  importance  that  could  be  incorporated 
in  the  model  is  as  fellows: 

(1)  The  capability  of  modeling  suspended  material  as 
well  as  that  in  the  bed  load  should  be  examined 
(cf.  MacDonald16). 

12)  Effects  of  using  a specified  distribution  of  sediment 
sizes  rather  than  assuming  a descriptive  mean  grain 
size  should  be  evaluated.  In  addition,  an  attempt 
should  be  made  to  determine  the  feasibility  of  ex- 
tending the  range  of  sediment  sizes  to  which  the 
model  can  be  applied. 


(3)  Accuracy  might  he  improved  by  modification  ox'  the 
assumed  forms  of  the  velocity  distributions  in  both 
the  laminar  and  turbulent  boundary  layers  (cf. 
Sleath55). 

I1*)  A further  evaluation  of  the  applicability  of  the  bed- 
load  function  over  a vide  range  of  conditions  would 
be  desirable. 

(5)  The  importance  of  using  nonlinear  wave  theory  to  cal- 
culate the  increase  in  wave  height  due  to  shoaling 
should  be  evaluated  (cf.  Jamesf0,71). 


PART  VII:  COHCLUDIHG  REMARKS 

1*7.  The  model  as  formulated  is  limited  to  transport  seaward  of 
the  breaker  zone  and  has  a number  of  simplifying  assumptions  and  limita- 
tions. It  should  be  potentially  valuable  as  a state-of-the-art  engi- 
neering tool  to  assist  in  the  placement  of  beach  nourishment  material. 
However,  a comprehensive  field  verification  should  be  performed  prior  to 
general  application.  The  importance  of  this  verification  cannot  be  over- 
emphasized; a comprehensive  and  controlled  data  collection  program  could 
serve  to  evaluate,  not  only  the  applicability  of  this  model  but  future 
ones  as  well,  thereby  permitting  potential  future  cost  reductions  for 
this  type  of  study.  The  principal  problem  encountered  in  attempting  to 
evaluate  applicability  of  sediment  transport  models  is  the  lack  of  reli- 
able prototype  data  for  verification. 

1*8.  The  tabulation  on  page  1*5  shows  estimates  of  longshore  sedi- 
ment transport  at  Point  Pedernales,  California.  The  absolute  estimates 
of  the  longshore  transport  predicted  by  the  model  (for  depths  greater 
♦han  12  ft)  are  about  15  to  22  percent  less  than  those  of  Bowen  and 
Inman**0  and  the  ratios  of  the  upcoast/downcoast  transport  are  quite  sim- 
ilar (0.83  and  0.90).  This  implies  that  while  further  refinement  and 
calibration  may  be  desirable,  the  physics  of  the  model  is  fundamentally 
sound. 

1*9.  Figure  10  shows  the  gross  annual  sediment  transport  predicted 
by  the  model  as  a function  of  distance  from  shore  and  water  depth;  no 
significant  longshore  transport  occurs  in  depths  greater  than  50  ft. 

The  dashed  portion  of  the  curve  indicates  the  zone  where  wave  breaking 
occurred;  therefore,  all  or  part  of  this  region  was  the  surf  zone  depend- 
ing upon  the  existing  wave  condition.  Figure  11  shows  the  net  annual 
sediment  transport  predicted  by  the  model,  again  as  a function  of  dis- 
tance from  shore  and  water  depth.  Note  that  at  a depth  of  about  20  ft, 
there  is  no  net  upcoast/downcoast  movement  indicated  and  that  on  an 
annually  averaged  basi3,  material  deposited  at  depths  of  20  ft  or  less 
can  be  expected  to  move  downcoast,  while  material  deposited  at  depths  of 
20  to  50  ft  can  be  expected  to  move  upcoast. 
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Figure  10.  Gross  annual  longshore  transport  as  a function  of  distance 
from  shore.  Point  Pedernaies,  California 
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APPENDIX  A:  NOTATION 

Amplitude  of  horizontal  displacement  of  water  particle  at 
bottom  due  to  wave  action 

Empirically  derived  constants  where  best  fit  the  theoretical 
curve  to  the  experimental  data  of  Kalk.au is  and  Abou-Seida 

Mean  stress  on  bottom 

Average  bed-load  concentration 

Wave  group  velocity 

Mean  sediment  diameter;  2D^q  » thickness  of  bed  layer 

Wave  energy  density 

Velocity  amplitude  function 

Phase  shift  function 

Acceleration  due  to  gravity 

Water  depth  below  still-water  level 

Wave  height 

Rate  of  transport  of  immersed  weight  of  sand  past  section  of 
beach 

Wave  number 

Von  Kerman's  constant 

Dimensionless  factor  of  proportionality 

A random  variable  that  has  a normal  distribution  of  zero  mean 
and  a standard  deviation  of  one 

Coefficient  of  proportionality 

Integer  determined  by  the  magnitude  of  the  Reynolds  number 

Ilondimensional  constant;  coefficient  of  proportionality 

Horizontal  mixing  coefficient 

Oscillatory  bed-load  rate 

Net  bed-load  concentration 

Reynolds  number 

Slope 

Friction  coefficient 


Wave  period 
Parameter 

Horizontal  velocity 

Maximum  orbital  velocity  of  waves 


A1 


Uj  Velocity  within  the  boundary  layer 


U(z)  Depth-dependent  horizontal  velocity 

Uw  Horizontal  fluid  velocity  outside  boundary  layer 

v Fluid  velocity 

vq  Maximum  horizontal  bottom  velocity  outside  boundary  layer 

V,X  Nondiaensional  variables 

Longshore  velocity 

v_  Average  velocity  at  a distance  Dc_  from  theoretical  bed 

v Mean  unsteady  velocity  within  boundary  layer 

x Horizontal  distance  normal  to  shoreline 

y Horizontal  distance  parallel  to  shoreline 

z Vertical  distance  positive  downward 

a Angle  of  wave  crest  with  respect  to  shoreline;  function  of  x 

Angle  of  wave  crest  at  breaking  with  respect  to  shoreline 

0 Wave  breaking  coefficient 

0Q  Reciprocal  of  boundary  layer  thickness 

Y Unit  dry  weight  of  bed  material 

s 

5 Thickness  of  the  boundary  layer 

6*  Laminar  sublayer  thickness 

c Roughness  diameter 

nQ  Empirically  derived  constant  which  best  fits  the  theoretical  i 

curve  to  the  experimental  data  of  Kalkanis  and  Abou-Seida 

A Phase  angle 

li  Dynamic  viscosity  i 

I 

we  Coefficient  of  horizontal  eddy  viscosity  i 

v Kinematic  viscosity 

€ Hiding  factor 

E Function  of  x 

Pf  Density  of  fluid 

p Density  of  sediment 

8 

pu  Density  of  seawater 

Ty  Met  stress  in  y direction 

+ Angle 

* Bed-load  function 


♦ Flow  intensity  function 

u Wave  frequency 


A2 


